Doping-induced metallicity and coexistence of magnetic subsystems in K 2 Fe 4+:r Se5 
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Electronic structure calculations are used to analyze the electronic and magnetic properties in 
K2Fe4+a;Se5. Fe atoms can be divided into two distinct groups. The x=0 (parent) compound forms 
an insulating, collinear, local moment phase with high Neel temperature. We show that large bi- 
quadratic exchange coupling and exchange-elastic interactions stabilize the magnetic order. For x>0 
the additional Fe atoms fill vacancy sites. They form impurity bands for small x, which broaden as x 
increases. They determine the states at the Fermi level and may be viewed as a magnetic subsystem 
separate from the host. Spin fluctuations are prevalent because magnetic interactions between the 
'defect' and the 'parent' atoms are relatively weak, while chemical fluctuations are prevalent for low 
x. Fluctuations of either type leads to the formation of a weakly metallic state. The unusual coex- 
istence of the two magnetic subsystems offers a new perspective as to how superconductivity and 
strong antiferromagnetism can coexist. We argue that spin fluctuations of the impurity subsystem 
share common features with the Fe-pnictide superconductors. 
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The iron selenide superconductors K^Fe^^Ses dis- 
covered recently^, demonstrate superconducting proper- 
ties similar to their pnictide counterparts, even though 
they appear to have very different magnetic ordering 2 
(compare Figs. [TJi, [lj>) and electronic structure near the 
Fermi level Ep£. This has a bearing on the associa- 
tions made between particular characteristics found in 
the pnictides, whose properties are largely invariant from 
one material to another, to the superconductivity ob- 
served. Spin interactions are likely to play a key role 
in mediating superconductivity; yet their character is 
not yet well understood^. While biquadratic exchange 
coupling 5 makes possible a consistent description of mag- 
netism in pnictides which form the striped antiferromag- 
netic (AFM) ground state (Fig. [Tf)), it is not clear what 
interactions are responsible for the stabilization of the un- 
usual magnetic structure observed in selenides 2 (Fig.[T]a). 
In addition, while weak (and most likely itinerant) AFM 
(Tjv~30-40K) and weak superconductivity (T G ~1-10K) 
have been observed to coexistence in the pnictides^, the 
situation is altered in this material. Magnetic moments 
and Neel temperature are much larger, strongly favoring 
the local moment picture, and the coexistence of AFM 
with superconductivity is much more robust (T)v~640K 
and T C ~30K)2&I. 

In spite of these differences, we will explain how it 
comes about that K2Fe4 +a; Se5 is metallic, how these two 
families share some key features in common, and further 
how these common characteristics can potentially pro- 
vide a framework to explain why these quite distinct ma- 
terial classes show similar superconductivity, while hav- 
ing quite different kinds of magnetic interactions. If such 
an hypothesis is correct, this work elucidates what the 
essential elements are for superconductivity to appear. 

The crystalline and magnetic structure of the parent 
compound K 2 Fe4Ses have been studied in RefA Vacan- 
cies substitute for 20% of the Fe sites in the unusual 



FIG. 1: Magnetic configuration of K2Fe4+a;Se5 in the block 
Neel structure. Parallelepiped shows unit cell. Fe atoms in 
the xy plane are shown: dark (blue) and light (green) de- 
pict spin-up and spin-down Fe atoms, respectively. There is 
a slight distortion of the square (sites 1-4) , elongating (short- 
ening) the diagonal along x (y). Open circles depict 4d sites 
that are empty for x=0, but get populated as x increases. 
Relaxation significantly stabilizes the magnetic order. Right 
panel shows same atoms in the striped AFM configuration 
found in the Fe pnictides, with x=l. 



ordered superstructure shown in Fig. [TJa. The Fe are 
arranged in squares of four atoms in a block with the 
spins parallel. The blocks themselves are aligned an- 
tifcrromagnetically, with two layers of atoms depicted 
in the Figure. K2Fc4Se5 has been identified as a mag- 
netic semiconductor— iSjJLS, with the Fe local moment ob- 
served to be ^i^B- It has been argued that supercon- 
ductivity appears in this system as a result of vacancy 
formation 2 ^^. 

For the parent system K2Fe4+nSe5, our (LDA) calcula- 
tion predicts a magnetic semiconductor in the block Neel 
structure with a local moment M=2.9 \xb and a bandgap 
of 0.44 eV, confirming the experiment and the findings of 
a prior study- . This latter work also concluded that lat- 
tice relaxation is necessary to stabilize the block Neel 
magnetic structure. We use the lattice constants of that 
work; our LDA implementation is described in Ref. 10 . 
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FIG. 2: Total DOS per Fe atom of K 2 Fe4 +:c Se5 in the collinear 
block Neel structure, for a;=0.25, 0.5, 0.75 and at x=l where 
all Fe vacancies are filled. The parent compound (x=0) is 
drawn in light grey, with a gap of 0.44 eV and the VBM 
aligned to 0. The remaining DOS are approximately aligned 
to the parent compound; Ep is drawn as a vertical dot-dashed 
line for each panel. As Fe begins to populate the vacancy sites 
(top panel), three distinct, very narrow bands appear between 
0.25 and 0.4 eV; Ef lies between the first and second. An- 
other localized state appears at _E=0.08 eV, and the VBM 
of the parent compound, splits. As x increases from 0.25 to 
0.5, the midgap states (especially the state just below Ef) 
broaden but begin to overlap, and when 1=0.75 they merge 
with the host bands. 



We consider how populating vacancy sites with Fe atoms 
affects the electronic spectrum and metallicity of the sys- 
tem, constrained for now to be in a collinear magnetic 
configuration. We denote Fe atoms in the parent com- 
pound as 16i atoms, and those filling the vacancy sites as 
4d atoms, following customary nomenclature. Supercells 
of the parent structure were generated, and a subset of 
the vacancy sites populated. In each case the lattice was 
relaxed to its minimum-energy configuration, using the 
PBE functional. (PBE-relaxed bond lengths, are in bet- 
ter agreement with experiment LDA ones; the average 
Fe-Se bond length in K 2 Fe4+oSe5 is 2. 444 A, close to the 
reported value. But we use PBE only to relax the struc- 
ture; the LDA is preferred for magnetic interactions.) 
Fig. [5] compares the evolution of total DOS D(E), with 
x against the parent compound, x=0. For x—0.25, a pair 
of localized states form in the gap: the lower band is filled 
and the upper band is empty. Both states broaden as x 
increases to 0.5; still they are nearly separated so that 
the system is at best weakly metallic (keeping in mind 



the LDA tends to overestimate hybridization and band- 
widths). Magnetic moments range from 2.1-2.3/xs (Fe on 
4d sites) to 2.7-2.9^s (16i sites) at 2=0.25. Thus in the 
absence of any fluctuations the impurity band widens, in- 
ducing a transition from semiconducting (D(Ep)=0) to 
weakly metallic behavior for x>0.25. 
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FIG. 3: Site and m-resolved d partial DOS D(E) per Fe atom 
in K2Fe4.25Se5. Ef is shown as a vertical line: the energy zero 
is chosen to align with the VBM of Fv2Fe4Se5. D(E) above the 
zero line are the m-resolved DOS summed over the 32 Fe 16i 
sites, while DOS below the zero line are summed over the two 
Fe 4d sites. The defect state slightly above Ef is confined to 
the host sites, while he one just below Ef, of 3z 2 — I character, 
is centered on the 4d site, with tails extending to the host. 

In Fig.^D(E) is resolved onto m components of d par- 
tial waves at the 4d and 16i Fe sites. The occupied defect 
band just below Ep is centered on the Fe 4d 3z 2 — 1 or- 
bital, with tails penetrating into the host whose cumula- 
tive weight approximately matches the head. The empty 
band, on the other hand, is almost completely confined to 
the host sites, and is analogous to a "surface" resonance. 

It is apparent that chemical disorder of Fe 4d atoms 
will lead to strong fluctuations in the local D(Ep): in re- 
gions where locally x> 0.5 defect bands will overlap and 
the system will be locally metallic. If the vacancy occupa- 
tion is random the local site occupation (concentration) 
will follow a binomial distribution, which can be reason- 
ably approximated by a Gaussian distribution. Then the 
probability P of finding a composition fluctuation in a re- 
gion containing N ions with amplitude exceeding Ai m i„ 
can be written in terms of the standard deviation a as 



P(Ax r , 



; N) = erfc ( A 



(1) 



where a 2 = x(l-x)/N. The defect wave functions are 
sufficiently short ranged that a sphere containing about 
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8 vacancy sites are sufficient to stabilize the local DOS. 
Partitioning an a compound with x=0.25 average va- 
cancy occupation into overlapping spheres of 8 vacancy 
sites each, about 10% of the spheres would contain a 
local concentration of x>0.5. Thus from percolation the- 
ory we can expect sizable portions of metallic behavior 
even at x—0.25, owing merely to local fluctuations in the 
concentration of 4d Fe atoms. 

Perhaps more important, there will be strong spin fluc- 
tuations on the 4d Fe sites, as we describe below. The 16i 
and 4d Fe atoms contribute to the electronic structure 
in approximately independent ways. The former have 
large local moments, are strongly exchange-coupled as 
described below, and form an insulating block magnetic 
Neel structure. They form relatively wide d bands with 
a large DOS; however, their contribution to the DOS is 
shifted away from Ep. Addition of 4d atoms causes de- 
fect levels to appear in the bandgap, which weakly couple 
to the 16i atoms to form a narrow defect band or reso- 
nance, at Ep. Though local moments of these atoms are 
calculated to be 2.3 /is their magnetic coupling to the 
host is weak. The reason is as can be seen in Fig. [TJa, in 
a collinear configuration, spins at the 4d sites must align 
parallel to one neighboring block and antiparallel to the 
other. The exchange field from the two blocks nearly can- 
cel, so the net exchange at these sites is weak. Therefore 
large, low-energy spin fluctuations must result. Such fluc- 
tuations can strongly reduce or even kill the Fe 4d static 
magnetic moment, making the 4d subsystem paramag- 
netic. This subsystem is thus similar to the entire lattice 
in other superconductors such as LiFeAs and FeSe. In 
such materials band theory predicts stable and large mo- 
ment while experiment clearly shows it is not present. 
In two recent Mossbauer studie s 11 ! 12 'nonmagnetic' Fe 
atoms were observed to coexist with the magnetic ones, 
e.g. in Ref . HP. 15% of Fe atoms in a K2Fe4.1_0.4Se5 sample 
were reported to be nonmagnetic. 

There are a host of recent experiments 2 -^ demonstrat- 
ing that antiferromagnetism and superconductivity coex- 
ist. This has prompted considerable debate in the liter- 
ature as to whether the two effects are present in the 
same phase, or whether a separate phase coexists that is 
responsible for superconductivity Several works present 
experimental evidence for a second phase; see e.g. Refi 14 . 
It is suggested that one phase is insulating in the block 
magnetic Neel structure that carries the antiferromag- 
netism, and the other in a nonmagnetic phase similar to 
FeSe, that carries the superconductivity There appears 
to be contradictory experimental evidence for both the 
"one-phase" and "two-phase" scenarios. 

The present work cannot rule out either scenario. But 
our findings show how it it is possible that magnetism 
and superconductivity can coexist in a single phase. The 
4d Fe atoms generate a rather itinerant channel or sub- 
system of states at Ep largely decoupled from the 16i 
magnetic structure. Fluctuations on the 4d subsystem 
will be large and the system as a whole weakly metallic 
even for small a; as a consequence of spin (and chemical) 



TABLE I: Average values of first and second NN exchange 
couplings, Ji and J2, in meV, for K2Fe4Ses. Interblock cou- 
plings are indicated with a prime, while intrablock couplings 
are unprimed. Also shown are the Neel temperature, esti- 
mated in the mean-field and random-phase approximations. 





M 


Ji 


J2 


Ji 


Jk 


rpMFA 

J-N 


rrRPA 

1 N 


unrelaxed 


2.9 


-5.7 


17.0 


25.7 


6.7 


682 




relaxed 


2.9 


-9.7 


10.2 


27.3 


8.5 


944 


494 



disorder. Separation of the 4d and 16i magnetic sub- 
systems is clear, as we show below, but the electronic 
states are somewhat mixed and eigenstates at Ep have 
projections onto both 4d and 16i atoms (Fig. Thus 
superconducting pairing still can originate from electrons 
of both subsystems, even while the metallic state is es- 
tablished by the small fraction of 4d Fe atoms, in a spin 
(and chemically) disordered configuration. 

Next we study the exchange coupling by using a 
Green's function linear-response technique, as described 
in Ref. 15 . Average values for such linear response J are 
shown in Table 1 for both ideal and relaxed F_2Fe4Se5 
structures. Exchange interactions more distant than sec- 
ond neighbor are found to be small, and we present only 
NN and 2NN parameters. We distinguish between intra- 
block couplings between FM aligned neighbors (Ji, J 2) 
and interblock couplings between AFM aligned neighbors 
(J^Jj). Relaxation stabilizes the magnetic order, as can 
be seen by inspection of the individual J's or from the 
Neel temperature estimated in the mean field and RPA 
approximations (Table 1). 

J 1 and Jj are radically different. This calls into ques- 
tion the customary interpretation of these parameters 
in terms of the Heisenberg model. We can extend the 
Heisenberg hamiltonian to include biquadratic terms, 
and assume^ that 

H = ^ JijSi • Sj' 7 Jij — Jij IKijSi • Sj. (2) 

ij 

For small angles, where linear response is applicable, 
can be eliminated if the intrablock and interblock are 
permitted to be different. It is apparent from the reduced 
symmetry that J need not be J'; moreover for small an- 
gles the anisotropic Heisenberg and biquadratic models 
are not distinguishable. They can only be distinguished 
at large angles; thus we evaluate the biquadratic term 
explicitly through calculations of the total energy E in 
large- angle noncollinear configurations. We accomplish 
this in practice by rotating the orientation of four atoms 
surrounding to one vacancy in a unit cell relative to the 
four surrounding the other vacancy (Fig. [Da), by a se- 
ries of angles 9 ranging between and w. Because this 
particular rotation preserves all second-neighbor angles, 
only NN terms in Eq. ((2|) contribute. A Fourier trans- 
form of E(0) yields directly the sum of all pairwise terms 
in Eq. ©, i.e. J = ^ J ij} K = V /; A,, and terms 
higher order in cos(9). Terms beyond the biquadratic 
are found to be small, so only J and K need be con- 
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sidered. As a check, J + 2K calculated this way should 
match J = J-ij calculated by linear response. Indeed 
we find this to be the case: the two calculations agree to 
within 3%. A large biquadratic coupling Kij of positive 
sign, on the same order as Jy, is necessary to explain 
the anisotropy in J and J'. This system is best char- 
acterized by large positive biquadratic coupling, which 
initially affects AFM coupling between blocks and FM 
coupling inside them. 

Lattice relaxation (which depends on x) stabilizes FM 
coupling inside blocks, as can be seen by its effect on 
Ji, Table 1. FM coupling is further stabilized by lo- 
cal distortions originating from partial filling of the 4d 
Fe sites. Bilinear and biquadratic exchange interactions 
both change with relaxation, even while magnetic mo- 
ment amplitudes are nearly constant. This effect, which 
we term the "exchange-striction" effect analogous to the 
well-known magnetostriction, is in part responsible for 
stabilizing the observed magnetic ground state. 

The magnetic structure and electronic states at Ep in 
iron selenides appear to be very different from the Fe 
pnictides. Nevertheless they share in common a large 
biquadratic coupling, which in each case helps to stabi- 
lize the magnetic ground state structure. In K^Fe^^Ses 
this mechanism is further affected by strong 'exchange- 
striction'. We showed that two nearly independent sub- 
systems coexist: a strong AFM phase with large local 
moments and high Neel temperature, and a more itiner- 



ant phase with fluctuating moments. Further, the metal- 
lic state (a prerequisite for superconductivity) is a con- 
sequence of fluctuations. This creates a plausible sce- 
nario to explain the coexistence of superconductivity and 
strong AFM observed in a single phase. In addition, 
both pnictides and selenides share two key generic prop- 
erties: antiferromagnetic interactions create a pseudogap 
in both cases, and the DOS at Ep is relatively low. This 
is also a prerequisite for spin-mediated superconductiv- 
ity because a large DOS at Ep can be a major source of 
incoherent spin scattering, which is usually destructive 
to superconductivity originating from spin fluctuations. 
The large moments we find on the 4d subsystem is likely 
an artifact of density functional theory — a mean field 
approach which does not incorporate fluctuations, and is 
expected to be very similar to the well known overesti- 
mate of the static magnetic moment in many iron pnic- 
tides systems. Thus the appearance of such soft 'itiner- 
ant' magnetic elements with weakly magnetic or param- 
agnetic spin fluctuations can be considered as the only 
generic feature of both iron pnictides and chalcogenidcs, 
while detailed electronic structure near the Fermi level, 
shape of the Fermi surface, magnetic structures and sym- 
metry of superconducting order parameter seems not uni- 
versal and not generic. 

Work at the Ames Laboratory was supported by US 
DOE Basic Energy Sciences, Contract No. DE-AC02- 
07CH11358 
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